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Abstract 

This  paper  presents  a  mathematical  model  for  the  coupling  between 
moisture  diffusion  and  damage  in  fiber-reinforced,  polymeric  composites. 
In  these  materials,  moisture  was  observed  to  cause  damage  by  a  multitude  of 
minute  debondings  at  the  fiber/matrix  interfaces.  The  model  employs 
concepts  of  continuum  damage  theory  to  describe  those  debondings.  Formal 
evolutionary  expressions  are  derived  and  related  to  the  extent  of  damage, 
the  stress  field,  moisture  content  and  moisture  gradient.  The  effects  of 
damage  on  moisture  diffusion  and  on  reductions  in  moduli  are  also 
formulated. 

Qualitative  comparisons  v/ith  experimental  results  are  provided. 
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Introduction 


It  is  wall  known  that  in  many  materials  deformation  under  loads  is 
associated  with  the  formation  of  a  multidute  of  internal  flaws.  These 
flaws,  which  may  be  micro  voids,  micro  cracks  or  micro  crazes,  precede  the 
development  of  macro  cracks  which  cause  final  failure.  The  abovementioned 
flaws  can  be  caused  by  environmental  agents  such  as  moisture  and 
temperature,  in  addition  to  mechanical  loads. 

In  the  many  circumstances  where  the  micro  flaws  are  distributed  in  a 
statistically  homogeneous  manner  it  is  advantageous  to  represent  them  as 
internal  state  variables  and  employ  thermodynamic  considerations  to 
establisn  constitutive  relations  and  evolutionary  expressions  for  flaw 
growth  ^  ^  1  *'  ^2] .  Tnis  approach  is  employed  by  several  "continuum  damage" 
models  which  were  reviewed  recently  by  D.  Kra jcinovic .  Guided  by 
various  ptiysieal  and  mathematical  considerations,  the  internal  state 
variables  were  chosen  to  be  scalars,  vectors,  and  tensors  of  various  ranks. 
The  case  of  a  vector  valued  internal  state  variable  was  employed  by 
Talreja^'^1  to  model  damage  in  fiber  reinforced  composite  laminates  and 
relate  stiffness  reductions  to  external  loads.  More  recently,  a  revision 
in  the  interpretatin  of  "damage"  as  micro-crack  areas  lead  to  the  selection 
of  internal  state  variables  as  axial  vectors  (or,  equivalently,  as  skew 
symmetric  tensors).^!  This  choice  will  also  be  employed  in  the  present 
work. 

Since  the  present  investigation  aims  specifically  at  fibrous 
composites,  where  damage  forms  in  characteristic  patterns,  the  existence  of 

^Numbers  "in  brackets  indicate  reference  listed  at  the  end  of  this  piper. 


a  "representative  damaged  cell  is  assumed.  The  components  of  the  axial 
vector  which  represents  "damage"  are  then  defined  as  the  projections  of  the 
total  area  of  microcracks  contained  within  the  cell  on  its  "walls".  When 
those  projections  are  divided  by  the  respective  areas  of  the  cell's  walls 
the  measure  of  damage  is  non-dimensional.  The  representation  of  all  the 
microcracks  within  a  cell  by  a  single  axial  vector  certainly  obscures  the 
distinction  between  few  "large"  microcracks  ana  many  smaller  microcracks. 
However,  in  circumstances  when  damage  forms  in  consistent  patterns  sucn  a 
distinction  may  not  be  important  because  the  variability  in  microcrack 
sizes  is  likely  to  be  limited.  The  interactions  between  microcracks  witnin 
the  cell  will  certainly  depend  on  the  external  loads.  it  will  be  shown 
that  the  present  model  accounts  for  this  dependence  througn  stress- related 
damage- evolution  relations 

With  few  exceptions ^1  ,  most  existing  continuum  damage  formulations 
employ  linearization  in  the  damage  parameters  By  contrast,  the  present 
formulation  does  not  involve  series  expansions  in  the  damage  parameter  and 
is  not  limited  to  "small"  damage. 

In  the  presence  of  sharp  gradients  of  temperature  or  moisture  content, 
the  expansional  strains  may  be  highly  non-uniform  within  the  characteristic 
damaged  cell.  In  this  case  the  stresses  are  likely  to  vary  even  along  each 
of  the  individual  microcracks,  resulting  in  elevated  stress- intensity 
factors  at  the  microcrack  tipa  Within  the  context  of  a  continuum  damage 
theory  these  increases  in  stress-intensity  are  reflected  in  gradient- 
dependent  damadge  evoLution  relationa  Sucn  relations  are  also  considered 
in  the  present  paper. 

The  effects  of  moisture  in  polymeric  composites  were  investigated  over 

I  3  I 

more  than  a  decade,  h  comprehensive  review  which  appeared  recently1-' 
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listed  more  than  three  hundred  references  on  the  subject.  Damage  due  to 

moisture,  which  developed  as  debondinqs  at  the  fiber /matrix  interfaces,  was 

observed  by  several  investigators^  This  typical  form  of  moisture- 

induced  damage  was  attributed  to  the  presence  of  hygrophilic  chemical 

f9l 

agents  at  the  fibers'  surfaces  .  Since  the  epoxy  may  act  as  a  semi-permeable 
membrane,  the  high  concentrations  of  moisture  result  in  excessive  osmotic  pres¬ 
sures  at  the  interfaces,  leading  to  fiber/matrix  debondings. 

n8i 

In  another  study1  J  it  was  shown  that  epoxy  resins  absorbed  excessive 
contents  of  moisture  when  the  amount  of  curinq  agent  in  the  mixture  was  below 
stoichiometry’.  Since  it  is  plausible  to  assume  that  the  stoichiometry  of  the 
resin  would  change  in  the  vicinity  of  the  fiber  interfaces  it  is  conceivable 
that  the  interphase  regions  contain  excessive  levels  of  moisture,  which  cause 
interfacial  cracking. 

The  process  of  moisture  sorption  is  associated  with  a  thermodynamically 
"open"  system,  since  vapor  mass  is  being  added  to  the  material  volume  of  the 
composite.  This  process  will  be  accommodated  in  the  present  paper  by  con¬ 
sidering  a  hypothetical  vapor  reservoir  which  is  in  thermodynamic  equilibrium 
with  the  actual  vapor  contained  in  a  material  volume-element  of  the  composite. 

ri9i-r2ii 

This  approach  follows  the  ideas  employed  by  Biot  in  connection 

with  flow  through  porous  media.  If  should  be  pointed  out  that  in  spite  of 

the  similarity  between  Biot's  approach  and  the  present  formulation  the  two 

are  not  identical.  The  natural  internal  variable  in  Biot's  scheme  is  the 

rare  pressure,  while  in  the  present  work  it  is  more  suitable  to  employ 

moisture  content,  or  alternately  the  chemical  potential.  The  subtle  differences 

[  221 

between  the  two  formulations  were  pointed  out  by  Gurtm 
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2.  Basic  Equations 


Consider  a  solid  oody  B  occupying  a  material  volume  V  boundea  Dy  a 
surface  A.  Let  tne  solid,  of  mass  density  Ps,  absorb  vapor  tnrougn  its 
boundary  and  let  m  denote  the  vapor-mass  per  unit  volume  of  tne  solid. 
Also,  let  x  be  the  position  of  a  solid  mass  particle  in  tne  deformed 
configuration  that  corresponds  to  tna  place  X  in  the  undetormeu  state,  and 
let  .f,  cj_  and  _v  denote  fluxes  of  vapor-mass  and  of  heat,  and  the  velocity 
of  tne  solid  particles,  respectively. 

In  addition,  let  u  and  s  be  tne  internal  energy  and  entropy  densities 
of  tne  solid/vapor  mixture  per  unit  solid  mass,  and  let  c^.  and  T  denote 
the  components  of  tne  Cauchy  stress  due  to  mechanically  applied  ioaas,  and 
temperature,  respectively. 

A  proper  accounting  of  tne  state  of  the  soliu/vapor  mixture,  wnicn  is 
a  thermodynamically  open  system,  is  obtained  by  considering  each  element  in 

r\j 

thermodynamic  equilibriuai  with  a  reservoir  containing  vapor  at  pressure  p, 

't  a, 

density  p,  and  internal  energy  and  entropy  densities  u  and  s 

respectively  &0  )  ,  [23  j  ,  124  j  ^ 

Conservation  of  the  solid  and  vapor  masses  gives 


p  +  p  V  *v  =  O 

s  s  - 


m  =  -7  •  f 


Conservation  of  energy  over  B  reads 

air/r-sudv=  Jannjv idA-  / 


d  A  -  I  qi  ni  d  A 


-  /  p  4"idA-  /“fini 
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3.  Piste lbuted  Damage 

When  materials  possess  a  statistically  homogeneous  microstructure, 
their  mechanical  response  is  associated  with  the  creation  and  growth  of  a 
multitude  of  internal  flaws.  For  several  types  of  material  microstructure 
these  micro  flaws  develop  in  characteristic  patterns,  until  they  finally 
coalesce  to  form  a  localized,  dominant  crack  whose  growtn  leads  to  ultimata 
failure.  Some  characteristic  damage  patterns  are  shown  in  Figs.  1  and  2 
for  a  fibrous  composite  laminates  of  different  lay-ups.  125],  [26]. 
Patterned  damage  was  also  observed  in  concrete  and  in  ceramic  materials. 

In  the  abovementioned  circumstances  it  is  possible  to  relate  the 

distributed  flaws  to  a  characteristic  material  "cell"  and  express  the 

[27] 

damage  by  means  of  a  continuous,  internal  state  variable  .  Such  cells 
are  overlaid  on  the  damage  patterns  in  Figs.  1  and  2. 

For  damage  due  to  micro-cracking,  the  internal  state  variable  can  be 
selected  to  represent  the  projections  of  all  micro-crack  surfaces  on  the 
"walls"  of  the  characteristic  cell.  Since  areas  are  expressed  as  vector 
products  of  directed  line-segments,  the  present  choice  loads  to  a 
mathematical  representation  of  "damage"  as  a  skew-symmetric,  second  rank 
tensor  d^jj.  The  quantity  d^jj  may  be  viewed  as  non-dimensional,  since 
it  can  be  formed  by  dividing  all  projected  micro-crack  areas  through  the 
respective  areas  of  the  cell  walls. 

In  the  presence  of  hygrothermal  effects,  diffusion  and  damage 
phenomena  are  Likely  to  depend  on  gradients  of  moisture  content  and 

of  the  temperature  ■T/  'X1.  Upon  consideration  of  the  abovementioned 
characteristic  cell  it  is  possible  to  relate  the  latter  dependencies  to 
non-dimensional  gradients  m/ •  1  x  and  •  T/  ■ r-  ,  where  fq  =  x^/tq  (no  sum  on 
l) ,  with  Li  beinj  the  lenjtns  of  the  cell  sides. 


1:  A  radiograph  showing  the  pattern 

matrix  cracks  in  a  [(),  90, A5]s,  Gr/Kp 
linate.  (After  Highsmith  et.al.,  Ref 
]).  The  "Character  i  st  ic-Uamage  (.'ell" 
super  imposed . 
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Finally,  it  should  oe  nocod  tnat  each  micro-crack  is  contained  within 
two  equal  and  opposite  surfaces.  Consequently,  tne  constitutive 
formulation  which  employs  d^jj  -as  an  internal  state  variable  should  remain 
insensitive  to  the  sign  of  d^jj. 


4.  Elastic  Response  with  Distributed  Damage. 

Consider  the  response  of  elastic  materials  with  distributed  internal 
damage.  In  addition,  let  the  material  be  expoed  to  thermal  effects  and 
absorb  moisture  from  the  ambient  environment. 

In  those  circumstances  the  list  of  internal  state  variables  contains 
the  deformation  gradients  F1j  =  '',Xj/gXj,  "damage"  dj^jj,  moisture  rn, 
temperature  T  and  the  gradients  Zj_  =  and  gx  =  gT/'.xi  of  the  chemical 
potential  and  of  T.  As  noted  earlier,  both  gradients  and  d^jj  may  be 
viewed  as  non-dimensional,  while  FiJ  is  obviously  dimensionless.  In 
addition,  T  and  can  be  non-d  imensional  i  zed  as  well  by  dividing  their 
actual  values  through  some  reference  levels. 

Considerations  of  frame  indifference  and  employment  of  the  reduced 
entropy  inquality  ^3)  ,give 


(ekl'  d[pq],  T) 


(3) 


mi 


’  'KL 


(9a) 

(9b) 

(9c) 


V^Y'iri 


where  .  =  u  -Ts  is  the  Helmholz  free  energy,  =^-(F^K  F^  -  6KL)  are  the 
Lagrangian  strain  components,  SKL  are  the  components  of  the  symmetric 
Kirchhoff  stress  ^9]  ^  are  the  components  of  the  damage  variable 

d^jj  referred  to  the  undeformed  conf iguration.  Since  both  d^jj  and  Djpgj 
are  areas  they  are  related  by 

d[ij]  =  J  xL,k  ekji  eLPQ  D [PQ]  (n> 

with  J  =  det  "x^/5Xk. 

In  addition,  in  eqn.  (10),  r[i^]  is  the  "affinity"  to  the  rate  of 
damage  growtn  '^>  I  i  j  J  =  d[ig]  '  nameay 


rr^i  "  f 


s  -d,  .  . 


Lij] 


Furthermore,  we  obtain  the  following  forms  for  the  fluxes 


(EKL' 

D[PQ] ' 

GC' 

9  ^ 9 

(EKL' 

D [PQ]  ' 

Gc  r 

Zg/  ^  f  nn) 

* 

v  [Id) 

(EKL' 

D  [  PQ) 

'  GC'  ZB' 

and 

'  IIJ] 

are 

components 

(13a) 

(13b) 

(13c) 


of  damage  growth-rate  in  the  reference  coordinates  Xy  They  are  related 
to  qx  and  fx  through  QA  =  XA  1  qL  Fy  =  XA  1  fx  and  !  jjjj  is  expressible 
in  terms  of  ‘  ^-jj  in  the  same  manner  as  given  in  eqn.  (11)  for  Djjjj  and 
d^jj.  In  addition  Gc  and  Xy  are  gradients  referred  to  the  undeformed 
configuration,  namely  g  =  x  i ,  c  'J  i  and  A\i  =  xi,b  zi 

To  simplify  the  subsequent  formulation  we  shall  restrict  ourselves  to 
isothermal  conditions.  in  this  case  g^  =  0,  whereby 

FA  =  FA  (CKL'  D l PQ] '  ZB '  ^  To> 


.-v  -a  • v 


(Ida) 
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and 

/N 

V[IJ]  =  '[IJ]  (EKL'  °[PQ]'  ZB'  To) 
In  addition,  the  inequality  (10)  reduces  to 


U,  .  f  .  +  r 
1  l 


[ijj  *[ij]  - 


(14b) 


(15) 


5.  Fiber  Reinforced  Materials.  Transvese  Isotropy. 

Consider  uni-directionally  reinforced  fibrous  materials.  Such 
substances  are  transversely  isotropic  about,  say,  the  x^  axis  and,  in  te 
absence  of  any  rignt-hanaed  or  left-handed  internal  structure,  possess  also 

•ft 

reflective  symmetries  in  the  xp  and  x^  axes  . 

To  derive  the  detaiLed  dependencies  of  v,  FA  and  :[jj]  on  EKL, 
and  ZQ  it  is  necessary  to  form  all  the  transversely  isotropic  invariants 
among  these  variables 

The  complete  list  is  given  in  the  Appendix,  where  A^j,  w^j  and 
denote  a  symmetric  second  rank  tensor,  a  skew-symmetric  second  rank  tensor, 
and  a  vector,  respectively. 

Note  that  the  Appendix  lists  33  invariants.  However,  the  thirteen 


invariants  I6,  I1L,  I L5,  I]_g,  I20'  I21'  l23’  *24'  x25'  l21>  x29'  :30'  x32 
are  odd  in  and  therefore  inadmissible  to  represent  damage  that  must 

remain  insensitive  to  the  sign  of  d  [  i  j  ]  -  This  leaves  20  invariants  for 
expressing  and  1 [ 1 j j ,  as  explained  in  the  sequel. 

In  view  of  eqn.  (8),  the  free  energy  .  depends  only  on  EKL  and  D[pg]« 
Consequently,  the  scalar  .  depends  only  on  the  ten  invariants  1-^  -  1^,  Ig, 

ho'  1 12'  119  dnd  l22’ 


These  combined  synmetries  are  denoted  by  the  class  T-4  in  reference  [30]. 


12 


Expressions  for  FA  are  obtained  by  considering  the  integrity  basis  tor 
one  symmetric  and  one  skew-symmetric  second  rank  tensor  and  two  vectors 
under  the  symmetry  T-4,  subsequently  retaining  only  those  terms  which  are 
linear  in  the  second  vector  The  results  of  this  procedure  are 
Listed  in  the  Appendix,  where  the  componets  of  the  vector  valued  function 


U:  U]_,  u2,  u3  ace  related  to  those  of  A^j,  W  ^ j  and  V y  The  thirty  two 
terms  -  P2g,  H1  "  t‘ie  exPcossions  for  (ui»  u2)  and  u3  3 re 
functions  of  the  20  invariants  1^  -  Ig,  l-j  -  I^q,  I^2  -  1 14 •  1 lyj >  1 19 r 
■^22'  *26'  log'  1 31  ond  1 3 3  formed  among  A^j,  Wjj  and  Vj. 

In  view  of  the  fact  that  the  flux  components  must  remain 
insensitive  to  the  sign  of  the  damage  variable  d^^j  it  is  necessary 
discard  all  terms  odd  inVJ^j,  hence  P3  =  P5  =  Pg  =  Pg  =  P^3  =  P33  =  P34  = 
p15  =  p17  =  P19  =  p20  =  p22  =  H3  =  H5  =  Hg  =  Hg  =  0.  This  leaves  only 
sixteen  terms  out  of  the  original  lists  of  32  terms  -  P2g,  Hj_  -  Hg. 


Expressions  for  are  obtained  by  considering  the  integrity  basis 

for  one  symmetric  second  rank  tensor,  one  vector  and  two  skew-symmetric 


second  rank  tensors  under  the  symmetry  T-4,  then  retaining  only  those  terms 


which  are  linear  in  the  second  skew-symmetric  tensor.  In  this  manner  a 


list  of  terms  is  generated,  which  contains  products  of  components  of  A^j, 
,4  and  with  one  component  of  the  second  skew-symmetric  tensor,  say 
At  tnis  stage  generate  a  second  list  of  terms  by  transposing  the 
indices  m  and  n,  namely  Dy  forming  scalar  components  of  transposed  terms 


among  A^j,  W  and  M x  that  correspond  to  Y^nTlj. 

The  desired  skew-symmetric  tensor  valued  functions  are  then  obtained 


by  suotr icting  the  factors  chat  multiply  Y^nn]  from  those  which  multiply 
^  t  run]  ’ 


The  results  of  this  procedure  ere  given  in  the  Appendix,  where  '’.[31]/ 
M32]  and  '{12]  ara  related  to  components  of  A^j,  W  ^  and  Vj.  The  eighteen 
terms  r^  -  r^,  h^  -  h^  in  the  Appendix  are  functions  of  the  same  twenty 

invariants  Ix  -  I5,  In  -  I10,  I12  -  I]_4»  Iig»  If7»  X19»  x22'  I26f  128t  :31 
and  I03  that  enter  P| ,  93,  etc.  Note  that  the  process  of  transposition 
from  Y[mn]  to  Y^nm]  and  subsequent  subtraction  automatically  eliminates  all 
terms  odd  in  /J  —  from  the  list  for  hence  no  further  reduction  is 
necessary  to  account  for  sign  insensitivity  to  the  damage  parameter  d^jj. 


6.  Infinitesimal  Deformations.  Strain  Formulation. 

Consider  infinitesimal  deformations.  In  this  case  Ejj  "•  ijr  skl'"  ~kl' 

Frkfi'  D [MN]  ^d[mn]  and  rs  "  >'so  (constant). 

Expanding  the  free  energy  in  powers  of  truncating  after  the 

second  power,  we  get 

1  so  =  :o  +  k  1  *33  +  ‘  2  (  "11+  c22)  +  ?3  C°(  :1l~f;22)  +  4d  [  3 1 J  d  [  32]  5  12^ 

2  2 

+  '  4d  [12]  (~31d  [32]  "  t32d[31])  +  Yle33  +  Y2(r~ll  +€22} 

,2 

+  '3  L’J(':11  -  22)  +  4d  [  31] d  1 32]  '12  -* 

2  2 

+  ”4  d  [  12]  (c31d[32]  _  r32d[31])  +  Y5c33(f:ll  +  ':22) 

+  ’6'  33  (  11  _r22)  +  4d{31]d[32]  C12^ 

+  '7  33d  [  12  J  ^ :  31d  [32]  “  >:32d[31]) 

+  ;3(:11  +  *  225  ^D(lll  “  e22^  f  4d[31]d]32]r12] 

+  'y(-  11  +  v'22)  d  1 12]  31d  [32]  “  f;32d[31]) 


g H* V ■  V «  V  tf  ■' If « V." w  l' «  IT  ',«n  '."J 1 ,  u  u  ',-v "V ■  »-vv 


B 
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'in  [Dl  = 


'10  LUV ’ll  “  l22)  +  4d[31]d[ 

2  2  2  2 


[ 31] a [ 32J ‘ 12^  d [ 12]  ^  ’  31d [ 32]  ~  f  32d [ 31]  ^ 


+  I 


11  ^Lli  ”  "22^  +  4t  i?  ^  f  'r;  ( 


12  J  +  '12  '31  +  "32) 


:  13d  [  12  ]  E^'ll  "  c22)  ('31d[32]  +  'V32d[31])  “  2’ 12 ( c3 ld [ 31 J  _ 


t'32d[32])  1 


(16) 


where  '• 'Qf  ^  (i  =  1,..#4),  r  L  (1  =  1,...13)  are  functions  of  m,  c3  ^  3 j 
2  2  2  2 

+  d  L 32 ]  ’  d [ 12] '  and  To-  Aiso  D  =  d [31]  ~  d  [32] * 


Stress-strain  relations  are  obtainable  from  o  .  =  p 

1 J  i)0 


,  wnare  £ 


11 


il 


is  considered  independent  of  t— .  This  requires  to  reconsider  the 
expansion  (16),  which  is  "biased"  in  favor  of  i:3i  and  12  and  r2Piace 

those  shear  strains  by  vj  (•,..+ c  ..) .  Upon  performing  this  modification, 

2  ig  jl' 

and  tnen  employing  the  "truncated"  notation,  withe Cy,  ’22  ^  °2' 


33'  3'  23 


31 


12  '*  '6 


and 


11 


22 


33  * '  3' 


2  23  “’4'  2--31*:  5 ,  and  2-;-L2'>-6'  we  °btai-n 


'p°  pq~q 


where  ^*pq  *“qp 


(17) 


In  (17)  ,  we  have 

C10  =  '  2+ '  3d'  C20  =  *2“ ''3°'  C30  =  Sl' 

C40  =  -2'4d[l2]d[3l]'  C50  =  2‘';4d[12]d[32] 

C60  =  2 '3a [ 3 1] d [ 32] 
and 

C11  =  2(  '2+,3u2+:3D+"'ll)  '  c12  =  2 ( *2-  1  3°2—  f  1 1 ) 


'13 


=  ; 5+  '6Df  C14  =2(-,9-  '10  D  +  '13^d  [  12]d  [31]' 


C15  “  2  (  ,9+V10d+'  13)dU2]d[32]  ' 
C16  =  2(2:3D+'3)d[31]d[32] ' 


-  *  -  *  -  *  -  •-  -  *  -  •-  -  —  a.  ~  -  -  a  —  j*  -  i  .  a  .1  -  *  -  a  a.  -1  «.  *<■  -  -  ►-  ~  w  -i  -  »  -  a  -  a  - 


Si 

V-is3 


C22  =  2(':2+",3D“-  '3D+  'll5  '  C23  "  V5~  '6°' 

_  1 

C24  "  2  ("  :9+  '10D~  'I3)d[12]d[31]  ' 

C25  =  2  ( V  i10D-  ,13)d[12]d[32]  ' 

C26  =  2  (“2  '3D+ !8) d  [31] d  [ 32]  '  c33=2tl'  C34  =  ~£' d [  L2 J d [  31]  ' 


'35  “  2  a [12] a [32]  '  ^36 

2,2.1 


'  2 '6  d [ 31] d [ 32] ' 


c44  -  : 4d [12]d  [31]  +  I  :  12' 


C 


45 


'’4d[l2]d[3l]d[32] 


-  —  1 1  ,~id  r ->  i  i  d  r  o')  i  d  m ->  i  +  d  f  -i  ,  d 


1 10u  [ 31] u  [  32]  [  12] 


2  2  ;  1 2 

u55  “  >4d[12]d[32]  +  -5-  '  C5G 
2  2 

c66  ~  2(4,3dL31]d[32]+*ll) 


[  12] a [ 32] ' 
( ! 10d  L  32] 


_  '13  , 


2  ' d [ 12] d [ 31] 


Note  that  when  all  d^j  vanish  tne  stress-strain  relations  reduce  to 

the  familiar  expressions  for  transverse  isotropy  about  the  x^  axis. 

However,  in  the  presence  of  damage  the  stiffness  matrix  Cp^  contains  all  21 

components,  all  of  which  depend  on  djjjj.  Obviously,  aLl  stiffness 

components  may  depend  also  on  m  and  Tq. 

For  1  we  neglect  all  terms  that  involve  the  symmetric  second  rank 

tensor  in  the  Appendix  and  obtain  the  following  expressions  for  the  vector 

that  represents  the  moisture  flux  _f: 

2  2 

fl  =  Plzl  +  P6l (d[31]_d[32] > Z1  +  2d[31]d[32] z2] 

2  2 

t2  =  Plz2  +  P61_(d[31]“dt32])z2+  2d  [  31] d  [32]  zl]  (18) 

f 3  =  H1z3  +  fI7dll2]  rd[31]  z2"d[s2]  zl> 

In  (18)  P|_ ,  Pg,  Hj_  and  H  -7  ate  functions  of  the  twenty  invariants  I  ^  - 

1 5 •  I7"I10f  I12_I14'  1 16 '  1 17 '  L19^  i22'  X26'  L28'  X31  jnd  *33  Eorn:ied  fro:T1 
the  components  of  the  symmetric  tensor  the  skew-symmetric  tensor  d^j 

and  the  vector  Zj.  Obviously  P^,  Pg,  and  H7  nviy  depend  also  on  the  m  and  on 


he  (constant)  temp-:  r  a  cure  Tri. 


Note  that  tne  expansion  of  .  in  power 


series  of  r  —  does  NOT  imply  that  an  analogous  expansion  must  exist  tor  P^, 
pb'  H1  and  h7'  Tde  Pcesent  formulation  therefore  retains  the  option  to 
consider  non-linear  coupling  between  mechanical  fields  and  moisture 
flux  t33l. 

Finally,  employing  similar  arguments,  we  obtain  the  following 
expressions  for  the  damage  growth  rates  =  d^jj: 

’l  12)  =  hld[12]  +  n4(d[31]z2  "  d[32Jzl)z3 

22 

31]  -  cld [31]  +  r4t(zl~z2)d[31]  +  z2(zld[32]  “  z3d[12])) 

+  riid[12]d[31j  zLz2  +  r  L3  (d  [  i2]  z3  -  d[32]zl)z2  (19) 

2  2 

.'[32]  =  c1d[32]  +  r4 1- (Zi-^2> d  [32]  +  zl(z2d[31]  "  z3d[21]) 

“  rlld [ 12]d [ 32] zlz2  “  r13(d[12]z3  _  d[13]z2)zl 

In  (19),  the  terms  h1,  h4,  rj_,  r4,  ru  and  rL3  are  functions  of  the 
same  invariants  as  PL,  p6,  h3,  and  H7  above.  Note  that  the  terms  and  r 
correspond  to  "seLf  similar"  damage  growth,  while  the  remaining  functions 
hg,  tg f  tgg  and  ri3  are  associated  with  the  "tilting"  of  micro-damage  due 
to  gradients  of  the  chemical  potential. 


7.  Infinitesimal  Derormation.  Stress  Formulation. 

Define  the  Gibbs  free  energy  d[ij]'  m'  To^ 

1  so  '  ‘so  '  Jijf ij 


(20) 


'hen,  in  analogy  with  cans.  (9)  we  have 


17 


Tne  "reduced  entropy  inequality",  eqn.  (10),  remains  unchanged,  except 


that  now 


*"  t  1 1  )  ‘  50  'dr. 

UdJ 

Consider  a  characteristic  material  stress,  e. g.  a  failure  stress  a £, 
then  for  sufficiently  small  stresses  (such  that  j  ^  j/-  £<<-l)  we  can  expand 
;  in  powers  of  and  truncate  after  the  second  powers.  This  expansion 
has  the  same  form  as  eqn.  (16),  except  that  replaces  j  ,  expansional 
coefficients  (-B^)  replace  3^,  and  compliances  (-rq)  replace  the 
stiffnesses  •,■.  An  analogous  procedure  yields  linear  strain-stress 
relations  similar  to  those  given  in  eqn.  (17) 

Lp  spo  +  Dpqcq 

with  3po  and  Spq  the  same  as  Cpo  and  Cpq  except  tnat  Bl  and  v i  appear  in 
place  of  .q  and  >q.  In  addition,  expressions  f  13)  and  (19)  also  remain 
unchanged,  except  that  P^,  Pfj,  ,  H7,  tq,  hq,  r^,  r^,  r  ^  and  r ^  depend 
on  twenty  invariants  that  contain  cqj  : n  place  of  iqj. 

In  view  of  (22b) ,  tne  chemical  potential  is  given  by 


qq  _  _Jo  _  R1  / !33\  ^2  /  711  +  C22  \ 

•n  '-m  In  \  '  /  '-m  \  / 

"  ^Td[l2]  ’[32]  (~— )  "  f1[3l]  (- 


rT  d[l2] 


+  signer  order  terms  in  ( 


stresses  1  he  flu;-:-  -s 


.  are  sr.vn 

i 


1  SC, 


m  •(.  -.  2  +  cj  c  ) 

[31]  [32] 


1  *  [  3 1  ] 


"•[31  ] 


« :¥) 


■  m  -dr 


\  d[l2] 


[.12]... 

,V 

‘'i  J 
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In  view  of  (25),  z.  will  depend  on 


^  '(d[3l]  +  d [32] ' 


,  and  on 


;  m  '•  d 


;m',(d[3l]  +  d[32]} 


3 

(i  =  l,  ...  4)  as  well  as  on  terms  like - ^ —  .  In  view  of  the  decendence 

•TO  ’X  . 

l 

of  It,  P, ,  . .  .r,  ,  in  eons.  (19)  on  .  it  follows  that  dr.  . -i  rav  introduce  a 
1  6  13  13  L13J 

non-linear  stress  effect  on  50  .  However,  for  sufficiently  short  times  -  when 


°[ij] 


Terrain  relatively  small  -  it  is  plausible  to  expect  that  for  n  . <<  1 


z.  will  be  linear  in  a. 

1  13 
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-VS.  ’V  -'v  . 


\  Special  Sub  Case:  Unidirectional  Diffusion  Under  a  Constant 
Transverse  Load!  "  ~ 

Consider  a  unidirectional ly  reinforced  plate  of  thickness  h,  witr 


floors  parallel  to  the  x-,  axis,  subjected  to  a  consca 


int  stress 


22  o 


with  diffusion  in  the  direction,  as  sketched  in  Fig.  3. 

In  view  of  tne  observation  that  most  damage  due  to  moisture  occurs  at 
the  f  ioer-matr  ix  interfaces  *  1  9  ]  ,  [  16  ]  ,  [  17  ]  ,  [  34  ]  assu:ne  j  =  o. 
Furthermore,  assume  gradients  only  in  the  direction-  Then, 
by  (2o)  ,  z2  =  Zj  =  0. 

in  these  circumstances  eqn.  (24)  gives 

h.  -  ^io  ^i2  o  ( ^  —  i »2 , 3 ,6)  ,  while  - ^  =  0. 

In  addition,  eqns.  (13),  (19),  (25)  and  (26)  yield: 

damage  growth  rates: 

a  2 

•  rr)i  -  —  C  ^  ^  i  ci  r  ->  ->  i  (27<i) 


■  [32]  "  Llu[32j  ~  I4"la[32] 

2 

[311  "  r i°  l  3 1  j  +  r4zld[31] 
moisture  fluxes: 

£1  =  plzl  +  p6d  Z1 

f 2  =  2p6d[31]d[32]zl 
chemical  potential  and  its  gradient 


(27a) 

(27b) 

(28a) 
( 28b) 


-  ar  -  -m(^)  *  £*($  *0«W2 
1  %  -  R2  \  +  >1  — +0<W2 


■i)( VV  (29b> 


where 


■2b 

2  1 

fjih 

•21 

2  ' 

•m 

Vc) 

rri 

2, 

111  ’(d[3l]  +  d[32]) 


■m  *  ( d  r  -j  i  i  +  d 


[31]  [32] 


T")  (t) 

nr  '  1  ' 


•■m  •U1[3]]  -  «[32] 
’’Fee  .also  Fig.  8  beli""w. 


•  (  ) 


In  tne  present  circumstances  rL,  t.,  Pp  and  P;j  m  eps.  (27)  and  (28) 

(  .  2  2  2  2 


depend  on  the  following  six  terms:  j  '0»  c'l3lJ  +  ^[321'  zl'  o  o  zl' 


zl  Df  as  well  as  on  m.  The  quantities  vQ,  E^*  and  B>  in  (29a) 


depend  on  dpi]  +  ^[32]  3nt^  on  ^ 


In  (28) ,  (29)  and  the  above  D  =  d. 


ill 


d[22]  as  before. 


In  view  of  (28)  and  (29b)  tne  process  of  moisture  transport  involves  a 


moisture,  stress  and  damage  affected  diffusivity.  Perhaps  more 

2  2 

significantly,  the  terms  -(dp-p  +  ^[32]^//':<l  -D/ 'x-|  in  (29b)  indicate 


that  sorption  is  influenced  by  damage  gradients  which  "cnannel"  moisture  in 
the  direction  of  increasing  damage.  Eqns.  (29)  indicate  that  c  depends 
linearly  on  stress  and  that  such  linear  dependence  is  likely  to  occur  also 
tor  the  diffusivity,  at  least  for  early  stages  of  damage  development. 

mined  by 

(30) 

where  is  the  chemical  potential  of  the  ambient  vapor.  For  small 


The  boundary  condition  on  moisture  content  is  deter 

h 


(X1  = 


2'  t}  ’’'A 


concentration  levels  it  is  plausible  to  assume  that  u  is  linearly  related 
to  m,  whereby  (29a)  predicts  saturation  levels  which*  at  least  for  early 
stages  of  damage  growth,  depend  linearly  on  the  stress  i  . 

An  experimental  investigation  of  stress-assisted  diffusion  in  AS4/3502 
g r apn i te/opoxy  coupons  was  con -Laded  recently  ^5].  unidir nationally 


reinforced  specimens  were  exposed  to  a  constant  relative  humidity  of  97a,, 
at  a  temp.ruture  of  40  C,  and  loaded  transvesely  to  tnv  fiber  directions  at 


0 , ,  15,,  301,  and  45.,  of  the  ultimata  stress  (where  'f  7500  ps  i  51.7 
MPa).  Total  moisture  weight- joins  were  recorded  p.-r  ioiical  ly  in  several 
tepl  irate  specimens  .and  resuits  for  tne  average  values  are  shown  in  Fig.  4 
below.  Note  the  "sigmoidal"  snnpo  of  all  absorption  curv  s,  v:h i  :h  differs 


qu  ilit  iti/ely  from  predictions  of  classical  diffusion  and  mdicit-s  »  non- 


vels  during  absorption 


linear,  concentr  ition  dependent  transport  process. 

Tna  dependence  of  the  maximal  moisture  content  end  of  tne  diffusivity 
on  stress  ire  snown  in  tigs.  5  and  6.  It  can  be  seen  that  on  approximately 
linear  reLationsnip  exists  between  stress  and  both  of  the  above  quantities, 
is  inferred  by  the  present  model. 

Moisture  weight- losses  were  measured  during  desorption  at  all  the 
joovement ioned  stress  levels.  These  measurements  were  performed  after 
removing  all  test  coupons  from  the  humid  chambers  into  a  dry  environment  at 
Os  relative  humidity.  The  resulting  weight  losses  are  plotted  vs.rfT in 
Fig.  7,  wnere  the  waignt-gain  data  are  superimposed  for  each  stress  level 
for  purpose  of  comparison.  Note  the  substantial  hysteresis  loops,  which 
can  be  attributed  either  to  the  concentration  dependence  of  the  transport 
process  or  to  the  growth  of  damage,  or  to  both. 

9 .  Moisture  Induced  Damage  In  The  Absence  Of  External  Stress. 

Consider  an  unstressed  unidirectional iy  reinforced  plate,  of  thickness 

h  as  before  with  ill  fioers  parallel  to  the  axis  and  moisture  diffusion 

in  the  Xj_  direction.  In  tnis  case  eqns.  (27)  and  (28)  remain  unchanged, 

2  2  2 

except  tout  r^,  r^,  and  depend  only  on  d  f  31  ]  +  ^[32]  '  Z1  anc^  rTU 
Eqns.  (29)  reduce  to 


3.00 


Fiqure  5.  Maximum  moisture  con 
stress  for  AS4/3502 


D  (10  7mm2/sec) 


0.2 


o.i 


Average  experimental 
values 

Approximate  relation 


0  0.15  0.30  0.45 

CTo/CTf 


I)  vii.:.  stress. 
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In  view  o£  eqns.  (27)  and  (31b)  it  is  cleat  that  damage  growth  rata 
depends  on  the  extent  of  existing  damage  and  that  this  rate  may  also 
depend  on  moisture  content  and  the  magnitude  of  the  moisture  gradient.  A 
dependence  of  vS[ij]  on  ^[ij]  wou^  lead  to  a  synergistic  effect  which 
accentuates  damage  localization  ("damage  breeds  upon  itself").  A 
dependence  of  upon  |7m/0xJ  wouLd  tend  to  localize  the  damage  in 

places  of  nighest  moisture  gradients,  namely  near  the  boundaries  x-^  =  +  h/2. 

The  typical  form  of  moisture  induced  damage  is  shown  in  Fig.  8.  Note 
that  "damage"  occurs  as  debondings  at  the  fiber/matrix  interfaces. 
Initially,  those  debondings  appear  as  isolated  inter  facial  crescents.  Upon 
repeated  absorption/desorption  cycles  tnese  crescents  grow,  until  they 
coalesce  to  create  highly  localized  damage  in  the  form  of  continuous 
cracks.  Typical  forms  of  such  cracks  are  shown  in  Figs.  9  and  10. 

The  growth  of  damage  can  be  inferred  also  from  weight-gain  and 
curvarture  measurements  in  anti-symmetric,  cross-ply  composite  plates^17'' 
136].  Due  to  the  anti-symmetry  of  the  lay-up,  those  plates  deform  into 
saddle  shaped  surfaces  upon  cool-down  from  the  elevated  cure  temperature, 
with  initial  curvatures  kx  =  -ky  =  k^.  Upon  subsequent  exposure  to 
moisture,  those  curvatures  vary  with  time,  whereby  k=k(t).  The  variation 
of  kx  -  k(t)  vs.  time  is  shown  in  Figs.  11  ana  12,  where  experimental 
results  are  compared  against  theoretical  predictions  of  linear  elasticity 
and  linear  viscoelasticity.  (In  those  figures  h  denotes  plate  thickness 
and  in  Fig.  11  tg  denotes  time  required  to  saturate  initially  dry  plat  is, 
prior  to  their  exposure  to  cyclic  ambient  humidities). 

The  variation  of  the  total  moisture  content  M  vs.  time  in  the  anti¬ 
symmetric  cross-ply  plates  is  shown  in  Figs.  13  and  14.  Weight-gain  data 
points  are  shown  in  comparison  with  predictions  of  classical  diffusion 


Time  ( days) 

Moisture  content  (in  7  weight 
duritir  eye]  je  exposure  to  07 
relative  tu iv  id  i  t  i  es  at  150'T, 
evele  interval  of  ')  da  vs. 


Inspection  of  Pips.  11-14  snows  increasing  departures  between  data  and 


theoretical  computations.  These  departures  are  most  likely  attributable  to 


Damage  Progression 


Shown 


i  us  1  ons 


A  continuum  damage  modei  was  developed  for  a  unidirectional ly- 
reintorced,  poiymeric-rosin  composite  chat  absorbs  moisture  from  a  humid 
amoionc  environment.  Damage  was  interpreted  as  the  total  cross-sect lona L 
area  of  microcrecKS  that  occur  within  a  characteristic  material  cell  prior 
to  the  formation  of  a  dominant  crude  The  total  microcracked  area  was  non 
d  i  mens  i  on  a  i  i  zed  through  division  by  the  respective  areas  of  the  cell's 
walls  and  was  represented  by  a  skew  symmetric,  second  rank,  censor  valued, 
interna!,  s  t  a  to  van  aol  e. 

M oistur-e  ingress  into  the  composite  was  treated  in  the  context  of  the 
tnermoiyn a.mcs  of  open  systems  and  coupled  moisture,  stress  and  damage 
relations  were  derived  from  fundamental  principles  of  thermodynamics  and 
continuum  mechanics.  These  relations  included  formal  expressions  for  the 
evolution  of  damage,  for  stress-and-damage-coupled  diffusion,  as  well  as 
for  u ana j e-dep._ nden c  material  compliances. 

It  was  shown  tnat  experimental  observations  of  moisture  induced 
damage,  and  of  moisture  absorption  and  desorption  in  the  presence  and 
uosenco  of  stress,  tended  to  verify  some  of  the  salient  aspects  featured  in 
the  continuum  damage  model  proposed  in  this  work. 

The  present  work  did  not  provide  expLicit  expressions  for  the 
evolution  of  d.amaje  or  any  of  the  other  formal  relations.  This  deficiency 
is  due  to  the  paucity  in  data  that  ire  available  at  the  present  time. 
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APPENDIX:  TRANSVERSE  ISOTROPY,  T-4  SYMMETRY  (about  X3  axis) 


Case  of  one  svnmetric  tensor  A. 

il 

one  anti-synmetric  tonsor  VX  , 

and  one  vector  V. 
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2.  VECTOR-VALUED  FUNCTIONS 
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